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A METHOD OF BIT AND POWER LOADING IN OFDM COMMUNICATION SYSTEMS 
WITH MODULATION AND CODING ADAPTATION 

Field of the invention 

5 The present invention refers to multicarrier communication systems, and more 

particularly it concerns a method of bit rate and power allocation to the carriers of one such 
system for transmission of each modulated symbol, and a system employing the method. 

The invention finds its preferred application in OFDM communication systems, and 
the following description will mainly refer to such preferred application. 

10 Background of the Invention 

Multicarrier communication systems are systems in which the band of a transmission 
channel is divided into a plurality of sub-carriers that are substantially independent and 
free of intersymbol interference (ISI). Symbols in the different sub-carriers are transmitted 
with suitable amplitude and/or phase modulation and suitable coding. 

15 One kind of multicarrier system is Orthogonal Frequency Division Multiplexing 

(OFDM), in which the frequency spacing between adjacent sub-carriers is chosen so as to 
make correlation between signals transmitted on the different sub-carriers null. An OFDM 
system allows using overlapping sub-carriers, so that efficient spectrum exploitation is 
achieved. 

20 Use of such systems has been proposed for high-speed data communication 

systems over both wired and wireless transmission media. For instance, in the wireless 
environment, a wireless local area network (W-LAN) using OFDM in its physical layer is 
being standardised by the European Telecommunications Standard Institute (ETSI). This 
network is known as HIPERLAN/2 and is disclosed in ETSI DTS/BRAN-0023003: 

25 "HIPERLAN Type 2 Technical Specification; Physical (PHY) Layer", November 1999. As 
known, in contrast with 3rd generation cellular system, which support high mobility and 
data rates up to some Mbit/s (e.g. 2 Mbit/s in UMTS), W-LANs offer much higher data 
rates (exceeding 50 Mbit/s) in a smaller range for stationary and portable devices with 
limited mobility. 

30 In an OFDM (or generally in a multicarrier) system, the problem arises of allocating 



bit rate and power to the different sub-carriers so as to maximise the overall bit rate for 
each symbol without exceeding a given power budget or, conversely, of distributing an 
assigned overall bit rate among the sub-carriers so as to minimise the required 
transmission power. 

5 Several proposals to this end are known in the art. Yet some of them are 

computationally efficient but sub-optimal, while other methods yield optimal allocation but 
are slow to converge and/or have a high computational cost. 

An exemplary optimal method is known as the Hughes-Hartog algorithm and is 
disclosed in EP 0 224 556 B. That algorithm starts by assigning no power and no bits to 

10 each sub-carrier. Then one bit or a fixed number of bits are given to the sub-carrier 
requiring the least amount of power to increase its data rate. The process is iterated until 
the power budget is exhausted. The Hughes-Hartog algorithm requires an intensive 
amount of sorting and converges very slowly in practical OFDM scenarios. 

Another optimal method, known as the Krongold, Ramachandran and Jones (KRJ) 

15 algorithm, is disclosed in US 6 400 773 B. This method exploits Lagrangian multiplicators 
to solve a constrained-maximum problem, where the function to maximise is the bit rate 
and the constraint is the budget power. The KRJ algorithm is cheaper to implement than 
the Hughes-Hartog algorithm, but it is actually optimal only if the allotted number of 
information bits is a convex function of the signal-to-noise ratio. Yet, this situation occurs 

20 usually when a single channel code and standard signal constellations (e.g. QAM, 
Quadrature Amplitude Modulation) are used. When different coding schemes and/or 
different modulation sizes are used, this may not be true and, under such conditions, the 
KRJ algorithm becomes sub-optimal. 

Object of the Invention 

25 Thus, the invention aims at providing a bit and power loading method that, while 

being computationally efficient, are optimal also when the coding and modulation lead to a 
non-convex bit-rate versus power curve, as it may happen when an adaptive coding and 
an adaptive modulation are used. 

Use of adaptive modulation in OFDM systems is well known. For instance EP 1 187 

30 413 A discloses a loading algorithm in which different modulation sizes can be used for 
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different blocks of sub-carriers. 

Summary of the Invention 
In a first aspect, the invention generally concerns a method of allocating 
transmission parameters (namely, bit rate and power) to individual carriers of a multicarrier 
5 communication system, preferably an OFDM system, so as to determine, for each symbol 
to be transmitted, an optimum overall value of a first of said transmission parameters while 
meeting predetermined constraints for an overall value of a second parameter. In case of a 
system providing for an adaptive coding and modulation, the method comprises the steps 
of: 

10 a) determining, for all carriers, a respective plurality of transmission modes characterised 
each by a specific combination of a modulation format and a coding scheme and 
consequently by a specific combination of values of said first and second parameters; 

b) determining, for all carriers, transition costs associated with a transition from one 
mode to a lower one characterised by lower values of said parameters, a transition 

15 cost being representative of a decrease in the first parameter as a function of the 

decrease in the second parameter; 

c) allotting to each carrier a respective first transmission mode associated with a 
maximum value of said transmission parameters; 

d) evaluating the overall value of said second parameter; 

20 e) if the overall value of the second parameter meets the constraints, using said 
respective first transmission mode for all earners; 

f) if the overall value of the second parameter does not meet the constraints, changing 
the mode allotted to one of the carriers to a lower mode, the carrier and the target 
mode selected for such change being those for which said change results in a 

25 minimum transition cost, 

g) evaluating the overall second parameter after the mode change at step f); and 

h) possibly iterating steps f) and g) until the second parameter meets the constraints. 

According to a second aspect of the invention there is also provided a multicarrier 
communication system comprising, at a transmitting end, means for coding and 
30 modulating a sequence of digital signals so as to obtain a sequence of symbols that are 
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each transmitted by using a plurality of suitably spaced carriers. The coding and 
modulating means are associated with a control unit, which adapts the coding and 
modulation laws and allocates bit rate and transmit power to the individual carriers by 
using the method of the invention. 
5 Brief description of the drawings 

The invention will be better understood from the following description of a preferred 
embodiment, given by way of non-limiting example, with reference to the accompanying 
drawings, in which: 

- Fig. 1 is a schematic block diagram of an OFDM communication system employing the 
10 invention; 

- Fig. 2 is a general flow chart of a first embodiment of the method of the invention; 

- Figs. 3A and 3B together are a more detailed flow chart of the method of the invention; 

- Figs. 4 and 5 are diagrams illustrating a simplified example the method of the invention; 
and 

15 - Fig. 6 is a general flow chart of a second embodiment of the method of the invention. 

Description of the preferred embodiments 
As shown in Fig. 1, an OFDM communication system employing the invention 
comprises a pair of identical transceivers 1 and 2 which exchange data through data lines 
3a (from transceiver 1 to transceiver 2) and 3b (from transceiver 2 to transceiver 1). As 

20 shown in greater detail for transceiver 1 , transmitter 4 includes: a coder 5 which receives a 
flow of information bits from a data source, not shown; a modulator 6 which receives the 
coded bits from coder 5 and provides symbols having phase and quadrature components 
(i. e. symbols of the kind l+jQ); a unit 7 computing an IFFT (Inverse Fast Fourier 
Transform) of the symbols generated by modulator 6 and transmitting the transformed 

25 symbols over data line 3a as a plurality of sub-carriers having the frequency spacing 
required by OFDM. 

According to the invention, both the coding law and the modulation law are adaptive 
and can be different for each sub-carrier. A control unit 8 determines the most suitable 
coding and modulation laws and, consequently, determines the manner in which the 
30 available bit rate and power are to be shared among the different sub-carriers. To this end, 
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control unit 8 receives from the remote transceiver 2, through a signalling line 9a, 
information about the channel conditions and, from upper layers of the system (arrow 11), 
information about the minimum bit rate and maximum power utilisable. Control unit 8 
performs the bit and power allocation algorithm according to the invention. 
5 A receiver 10 in transceiver 1 will receive OFDM symbols from transceiver 2 through 

data line 3b and will provide the remote transmitter, through a signalling line 9b, with 
channel information allowing a remote control unit to determine the most appropriate 
coding and modulation laws and performs the bit rate and power allocation for the signals 
to be transmitted over line 3b. 

10 The theoretical principles of the algorithm will now be disclosed. Without loss of 

generality, we assume that no inter-symbol interference (ISI) nor intercarrier interference 
(ICI) exists between the N sub-carriers of the OFDM system, so that the communication 
channel can be decomposed into an equivalent set of N additive white Gaussian noise 
(AWGN) parallel channels. 

15 Moreover, we assume that the number of information bits (modulation size) as well 

as the channel coding scheme and the allotted power can be adapted for each sub-carrier. 
We will thus use the term "transmission mode" to indicate the combination of a given 
channel coding scheme and modulation format for a sub-carrier. 

Let b[m] and b c [m] be the number of information bits and coded bits per OFDM 

20 symbol, respectively, on sub-carrier m, with m=l,...,7V . Hence, the code rate on sub- 
carrier m is R c [m] = b[m] / b c [m] . We observe that while b c [m] determines the 

modulation type by assigning the number of coded bits carried by each modulation 
symbol, R c [m] determines the coding scheme. Last, let P[m] be the transmit power on 
sub-carrier m and P = [P[V,P[2],...,P[N]] the vector of allocated powers. The overall 

25 number of information bits conveyed by an OFDM symbol is given by 
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and the total power is 

We indicate the set of available transmission modes as M = {0, 1...K} . Each mode is 

JV 

p= 5Z p H- (2) 

identified by a couple of numbers (bk,b Ct k ),k e M, where Jty e B = {6q = 0 f i>| ,...,/>*} and 
b ck eB c = {b c0 = 0,6 c1 ,...,b cK }, indicating the available number of information bits per 

5 sub-carrier use and the number of coded bits per sub-carrier use, respectively. Hence, the 
code rate of mode k e M is R c ^ =b^/ b c ,/c . In general, the first mode for k = 0 

corresponds to no transmission and we set R c q - 1 . Moreover, set B is ordered such that 

We denote with M = [M[1]... M[1]] an /V-uple of selected modes. 
10 If p e [m] is the bit error probability of sub-carrier m, we assume that a target error 

probability P e ,ref must be met in each sub-carrier, which probability is usually the same for 
all sub-carriers. I.e. it must be Pet'"] <Pe,re/> m = 1—N. 

Optimal bit and power allocation requires solving either of two maximisation 
problems, namely: 

15 1) Throughput maximisation. Given a fixed amount of power P to t , the objective is to find 
the A/-uples M and P that fulfil 

max b 

P = Puu < 3 > 
P e \m] < Pesef > Vm 

2) Margin maximisation. Given a fixed throughput bf 0t , the objective is to find the /V-uples 
M and P that fulfil 

ruin P 

M<= H v t P 

b = b tot 



As known, "margin" (or better "performance margin") is the amount of noise the 
system can tolerate while still operating under the bit error probability constraints: 
minimising power, as indicated by (4), actually corresponds to maximising margin. 

With regard to sub-carrier m, we denote with G[m] its gain and with <r* the noise 

5 power, which we assume to be the same for all sub-carriers; then we define the 
normalised signal-to-noise ratio for a transmit power equal to one, as 

Thus, for a transmit power P[m] , the signal-to-noise ratio is readily given by 

r(m) = P[m]f[m] (6) 
We will consider in greater detail the throughput maximisation problem. The 
10 considerations made for this problem can be repeated substantially unchanged for the 
margin maximisation problem, once the roles of the bit rate and power are exchanged. 

A. Throughput Maximisation Algorithm 
The constraint on the error probability is such that each transmission mode 
corresponds to an operating point in the b[m] -T[m] plane. We indicate with 
15 Tfo,k e M,(Tq =0) the signal-to-noise ratio required to transmit in an AWGN channel with 
a bit error probability P e ,re/> using mode k. We now introduce various metrics associated to 
the generic sub-carrier m. Firstly, for a given transmission mode k, the minimum power 
level which guarantees the target bit error probability is given by 

P fc fm]=J^T f *€M. (7) 
T[rn\ 

We also define the bit-per-power cost function corresponding to the transition from 
20 mode k to mode h: 

For a given mode k, we evaluate the mode x^fm] , lower than /c, which yields the 
lowest cost function. In other words, given mode k, h = x^[m] denotes the mode that 
yields the smallest information bit reduction ( - ) for a corresponding power 
reduction P[<[m] -P^[m] . I. e., in a plane (b, P), segment k - h is the lowest slope 
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{arg 



arg miu JF^fm] , kj^O 
x fc \m] = { h <* . (10) 

fc = 0 



segment among those connecting node k to any lower node h. 
Let 



{miu 



then 

Moreover, for each sub-carrier m, let M[m] represent the current allotted mode, X[m] 
5 the next optimum lower mode and F[m] the cost function associated to mode transition 
M[m] -> X[m] . 

The basic idea of the algorithm is as follows. Reference is made also to Fig. 2. First, 
we allocate the mode with the highest number of information bits to all sub-carriers, i.e. we 
set M[m] = K (step 20 in Fig. 2). This gives the maximum total number of bits b = Nbf < . 
10 Correspondingly, the initial power allocation is P[m] = Pfc/Vn/, for m = \...,N . For sub- 
carrier m, we also select the next optimum cost function F[m] = Fx[m] with corresponding 
mode X[m] = xk [m] . 

For a given power allocation, we then evaluate the total power P. If the power 
constraint in (3) is not met, we change the mode in the sub-carrier with the lowest cost 
15 function, i.e. set 

fh = arg nun Ffm) , (11) 



it is 



M[m) = X {rn] . (12) 
We then update metrics of sub-carrier m . 

Steps (1 1) - (14), which correspond to step 21 in Fig. 2, are repeated until P < P tot . 

It is easy to show that the "mode reduction" operation (12) is optimum in the sense 

P[m) = P M | m ,(m] t b[in] = b M|mj ( 13) 

F(mj = 7" M |m|N : Xfeft] = x M] , n \[™) - (14) 
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that, taking the total number of information bits and overall power after reduction as 
reference, we cannot transmit more information bits with the same amount of power or, 
equivalently, we cannot transmit the same number of information bits with a lower amount 
of power. In fact, the algorithm reduces the number of information bits by removing one 
5 mode at a time, discarding the modes with the "most expensive" bits, i.e. those bits that 
require the greatest amount of power to be transmitted. 

We observe that the algorithm works by "mode removal" rather than "bit removal". 

We note that this procedure is optimal for either throughput or margin maximisation. 
The only difference between the two problems is the constraint that must be fulfilled (a 
10 power constraint in throughput maximisation, a bit rate constraint in margin maximisation), 
which leads to a different stopping rule of the mode reduction operation. Taking this into 
account, the appended claim 1 and the claims depending thereon will generally refer to 
maximisation of a first transmission parameter while meting certain constraints for a 

second transmission parameter. 
15 In the particular case the number of information bits of mode k is a convex 

function of (convex case), i.e. the operating points in the - r> plane belong to a 
convex hull, metric (9) becomes 

r 



and relation (10) simplifies into 



(15) 

oo . k = 0 



■I: 



- 1 , k ^ 0 

arjbH ={ - ( 16 > 

k=0 



This means that the transitions are always between one mode to next lower one. 
20 We recall that usually b^ is a convex function of r> only when a single channel code 

is used and the transmission modes employ a standard signal constellation (e.g. QAM), 
whereas this may not be true when different coding schemes and/or different modulation 
sizes are used. 

At the end of the mode reduction procedure, if the total allotted power is below the 
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budget, i.e. P < P tot , and b < Nb k , then the choice (11) may be no longer optimal 
because, for some sub-carriers, there could be another transition that removes fewer bits 
and still meets the power constraint, although requiring more power than the previous 
selected mode. Hence, these transitions that keep the power below the budget P tot 

(hereafter named "cross transitions") do not satisfy (9). On the contrary, from the above 
discussion, all transitions that keep the power above the budget satisfy (9). As a 
consequence, at the end of the procedure, we have to compare all possible cross- 
transitions associated to the sub-carriers and choose the one that yields the greatest total 
number of bits b (step 22 in Fig. 2). 

Overall, the throughput maximisation algorithm has the following steps, where M opt is 
the optimum mode allocation, with corresponding number of bits b op t . Reference is made 

also to Figs. 3A and 3B. 

1) Initialise M = [K,...,K],P = tP k [V,...,P k [N]],b = Nb k ,P = z"^ P k [m] (Step 31, Fig. 
3A) 

IfP <P tof (Output Yes from Step 32, Fig. 3A) then M opt = M and b opt = b (Step 33, Fig. 
3A). Stop 

Else (P>Ptot) ' set b opt = 0 and initialise F[m]= F^m] and X[m] = x k [m] , for 
m = l,...,N (Step 34). 

2) If P> P tot (Output YES from Step 36): 

If F[m] = oo , for m = 1,..., N , (Output YES from Step 35), then stop. 

Else : 

• Reduce mode according to (11) - (14) (Step 37) and compute total number of 
information bits b and total power P. 

• Go to step 2), i. e. return to Step 35, fig. 3A 

Else (P <P t ot) (Output NO from Step 36, Fig. 3A) 

IfP- Ptot (Output NO from Step 40, Fig. 3B), then M op * = M and b opt = b. (Step 41, 
Fig. 3B). Stop 

Else (P < Ptot) (Output YES from Step 40, Fig. 3B) 

• ifb> bopt (Output YES from Step 42, Fig. 3B), then M opt = M and b op t = b. 
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• Back-off (Step 44, fig. 3B): undo the mode reduction in the last selected sub- 
carrier m . Let M[m] and X[m] be the new metrics (with corresponding 

P>Ptot)- 

If X[m] <M/m;-1 (i.e. there are still modes to be visited), select the next 
5 mode, not yet explored, with a higher number of information bits than any 

previously selected modes, and with the lowest cost function among the 
remaining modes. Set x[m] to this mode and F[m] = F M[m]X[m][™J 

Else ( X[m ] = M[m 7 - 1 ) (i.e. there are no more modes available), then set 
X[m] =M[m]- J \ ax\6F[m] = oo. 

10 •Goto step 2 (i. e. return to Step 35, fig. 3A.). 

Note that before the back-off operation, M opt is updated only if the last mode 

reduction yields a number of bits b higher than any previous allocation. 
Typically, the algorithm stops when F[m] = oc^ for m = 1,..., N . 

If desired, and as disclosed at step 45 in fig. 3B, the allocated powers can be 
15 normalised such that the total transmit power is equal to the budget P to t . 

In practically carrying out the invention, some measures allow reducing complexity. 

First, instead of pre-computing and storing cost functions (8), we can save 
computational effort and memory locations by calculating and storing in a matrix c (see 
example of Table I) the following normalised cost functions, that are independent of sub- 
20 carrier index m 

F'kh = 7T - * h . k = l,...>K . h=Q K-l t k>h. (17) 

In the convex case, for each mode k, we need to store only value F'^_^; 
otherwise, in general, we have to store values corresponding to h = 0,...,* — 1, which are 
sorted in increasing order. However, for each mode k, after sorting, we can discard 
transitions that lead to a mode with a lower number of information bits than any previous 
25 transition, because these transitions will never be visited in the back-off operation. Using 
(7), we can calculate the cost functions (8) from (17) by a simple multiplication 

jpfcftH = Ft* ?H . (is) 

Moreover, we note that X^[M] in (10) depends only on F' and not on f[M]. This 
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can be easily seen by substituting (18) into (10). 

Actually, along with matrix c, we can build up a matrix x (see example of Table II), 
which for each mode k stores the "arrival modes" indicated by the corresponding matrix c. 
In particular, in the convex case, only one value is stored per row, namely x[/c,1] = k-1, 
5 /c = 1,...,K. 

As an example, Figs. 4 and 5 show functions b k versus l~ k and b k versus P k (m), 
respectively, for a system of three sub-carriers, sorted in increasing order of f, and five 
modes. The transitions of interest for the method are shown in Fig. 4 and for one carrier in 
Fig. 5. Tables I and II show how to build the matrix of normalised cost functions and the 

10 matrix of arrival modes, respectively, and hence how to trace Figs. 4 and 5. 

The Figures clearly show transitions (the "cross-transitions" mentioned above) that 
have higher costs F' and yet entail a lower bit rate reduction than a lower cost transition. 
Such transitions are considered in the back-off operation. For example, from mode k = 4 
we can go to mode 0 or, with higher cost, to mode 3. A transition to mode 2 would also be 

15 theoretically possible, yet it would result in a higher cost and a lower number of information 
bits than transition to mode 2, so that it will never be selected in the back-off operation. 
Values corresponding to transitions that will never be visited by the algorithm have been 
omitted in Tables I and II. 

Note also that, in a situation like that depicted in Figs. 4, 5, the KRJ algorithm would 

20 never take into account modes 2 and 1 , which are out of the convex-hull approximation of 
the curves b k - r k and b k - P k (m) and therefore the allocation cannot be optimal. 

TABLE ! TABLE II 

MATR3X C OF NORMALIZED COST FUNCTIONS MATRIX X OF ARRIVAL MODES 



Mode 1 


n 0 








Mode I 
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Mode 2 










Mode 2 


0 






Mode 3 




r 31 






Mode 3 
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Mode 4 










Mode 4 
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Mode 5 










Mode 5 
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Additional saving in computation can be achieved by ordering the N sub-carriers 
such that the normalised signal-to-noise ratios \r[V~.r[N]\ are in increasing order. 

Moreover, we define a new metric vector, F' = [F[1]...F[N]], which contains for each sub- 
carrier m the normalised cost function f ' 
5 For each entry of matrix C, of value f'm,, we define a list containing all sub-carrier 

indexes whose normalised cost function is f in other words, this list points out all 
entries of F* whose value is F '*/,. Let L be the number of non-zero entries in matrix C. We 

introduce lists g (1) ...g (L) , where the superscript / (/ = 1...L) enumerates the entries in C in 
columnwise order. Hereafter the correspondence between f ' kh and / will be denoted by / = 
10 s(F k h ). For example, with regard to Table I, at a given step of the loading algorithm, g (1) 
will contain the sub-carrier indexes whose associated cost is F '10, g <2) those with 
associated cost F '20, and so on up to g (L ) which contains the sub-carrier indexes whose 
associated cost is f 32. We assume also that the elements of each list, g (,) , are sorted in 
increasing order. Hence, because the normalised signal-to-noise ratios T[m],m = 1,...,A/ , 

15 are in increasing order, the first entry, g (,) [1], denotes the sub-carrier with the smallest cost 
function F kh within list g (,) . 

We also introduce a vector, F"™, of L elements, containing the smallest cost 
functions F for each of the L lists. In the example, 

f-» = [^of^^liIl.^ofb^Wh-.-.^rbwii]]]. 

If a list is empty, the corresponding elements in F mln is set to 00 . 
20 We can now simplify the evaluation of (11). Let 

7= argminF rain [tj , (19) 

it is 

m = g (i) [l| . (20) 

A step-by-step description of the algorithm is given in Table III, where we assume 
that the normalised signal-to-noise ratios [r[y,...,T[N]] are sorted in increasing order. 

We also denote with pop(g 0) ) the operation of extracting the value from the top of list g 0) , 
25 pop(g Q) ,m) the operation of extracting the element of value m from g 0) , and push(g 0) ;m) the 
operation of inserting, in increasing order, the element of value m into g (i> . In Table III the 



vector J =[J[l],...,y[w]]contains indexes. 
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TABLE III 

THROUGHPUT MAXIMISATION ALGORITHM 

1) Initialize M = \K, . .. , K], P = (Pk[1], . . . , F*(A% b = JVbjr , P = 2Z™=i M- 

P < f\^t, then ^Upt = M and = P. Go to step J. 
£2w isetfi^ = 0,hiild matrixes C and A\ Initialize vectors J = [l f ...,l| r F* - [C|/C, . , Cl#C.j(A f l]]* 
X = [*[K,i[lJ], ... JMI, E <*> = [1, . .. , JV| and P«« = [f^I f (1 |,*o *>]. 

2) rf7&fefVoo 

ff P > P**t 

o Reduce made : 

i - arg mm F"**" [il 

m = g^[l| , pop(^) , F" pi -^liBFhP mi 



Pm = P 

i«u = j\nt\ 
fc[M|mJ t jfm|| , M[m| /o 



|*u , M[tfi) =0 

/=*(F^[m)) , push( E "\m) , tfg\ r > ^mitowFS? = F™fl aiid F^fl = F>|f [m| 



J*|M|m|, Jim]] , Mfmjj£0 
X|mJ = < 

^o, M[m\ - o 



o If 6 > fit,^, then - M, P„ p * = P and 6^ = 6. 
o Back-afT : 

M„c- = M|mj , M[m| = m^j , P[m| = PwjmjH f o[m| ■= b M |m] 

P = P^J i b = *oU 

;/g ,0 [l] = m <fce« F^ff) = F^j" 
pop( B (/) ,m) , j|m| = ;Mj+l 

J C|M[m|,i|m|| , -M„ ro ?M r fc J A-|M|m|,j|m|j . - M nt(r ^ I 

I = a(F'[m{) , push(s (n , m) , jTg fl> [l] = tn then F"" B [i] = F"[mJ f [mj 

(P = Ptot) : Mope = M and P^t = P Co to step 3. 
End While 

3) Normalize power : P, P » ftn] - P^, |ni] - m = I . . . . t A r . 



BEST AVAILABLE COPY 
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Point 3 is the power normalisation discussed above. 

Two possible simplifications can be applied to the algorithm to further reduce 
complexity: 

1 ) Convex case approximation : We can discard the modes that do not belong to the 
5 convex hull in the b-T plane. In this case only one transition from each mode to the 

next lower one is visited by the algorithm. This approximation may lead to significant 
computational saving with usually negligible loss in throughput 

2) No back-off : We can skip the back-off operation and terminate the algorithm as soon as 
P becomes lower than the budget P tot . Although the resulting allocation is no longer 

10 optimal, negligible throughput loss is expected as, usually, the back-off process may 
change the assigned mode in only a few sub-carriers. 

An evaluation of the complexity of the proposed algorithm shows that such 
complexity is substantially comparable with that of the KRJ algorithm, when applied in the 
conditions leading to an optimum solution also for the latter algorithm. We however point 

15 out still once that the present invention leads to an optimum solution also when the curve 
b(l") is not convex, whereas under such conditions the KRJ algorithm is only sub-optimum. 

8. Margin Maximisation 
We now consider the margin maximisation problem given by relation (4). Reference 
is made also to Fig. 6. As in the case of throughput maximisation, we start with allocating 

20 the maximum number of information bits and the maximum power to all sub-carriers. As 
clearly shown in the flow chart, the procedure is similar to that disclosed above for 
throughput maximisation, except that the number of information bits instead of the power is 
evaluated at each step to start or to go on with the mode reduction or with back off, 
respectively. In this case, reducing the rate in the m-th sub-carrier, according to (11M12), 

25 is optimum, provided that the corresponding total number of bits does not fall below the 
constraint b to t. If the target number of information bits b to t is not reached exactly through 
mode-reduction operations, the algorithm goes into back-off mode as in throughput 
maximisation and stops when b is the closest lower approximation of b to t. Indeed, it is to be 
appreciated that bit rate b is a discrete number, whereas power P considered in the 

30 throughput maximisation is a real (analogue) value. 
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The algorithm code is reported here 

1) Initialise M = [K,...,K], P= [P k [U...,P k [N]], b = Nb K , P = I^ =1 P k [mJ 
Ifb <to tof , then M op t = M and P opt = P. Stop 

E/se (/> > b tot ): set P opt =0 and initialise F[m] = F^m] and X[m] = x k [m], for 
m = 1,...,A/. 

2) If b> bfot 

If F[m] = oo , for m = 1,..., A/ , then Stop. 

E/se: 

• Reduce mode according to (11) - (14) (Step 37) and compute total number of 
information bits b and total power P. 

• Go to step 2. 

Else (b <b t od 

lfb = btot, then M op f = M and P op t = P. Stop 
E/se (to < to to ^ 

• ifP > P op t, then M op j = M and P op t = P- 

• Back-off (Step 44, fig. 3B): undo the mode reduction in the last selected sub- 
carrier m . Let M[m] and X[m] be the new metrics (with corresponding 

b>b tot y 

If X[m] < M[m] -1 (i.e. there are still modes to be visited), select the next 

mode, not yet explored, with a lower power than any previously selected modes, 
and with the lowest cost function among the remaining modes. Set x[m] to this 
mode and F[m] = f M[m]X[m] [™1 

Else (X[m] = M[m] -1 ) (i.e. there are no more modes available), then set 
X[m] = M[m]- J \ andF/7r?7 = oo. 

• Go to step 2. 

The skilled in the art will appreciate that the above disclosure is given by way of 
illustration only and has no limiting purpose, the scope of the invention being defined by 
the annexed claims. 



